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Frequency conversion (FC) is an enabling process in many quantum information protocols. Re¬ 
cently, it has been observed that upconversion efficiencies in single-photon, mode-selective FC are 
limited to around 80%. In this letter we argue that these limits can be understood as time-ordering 
corrections (TOCs) that modify the joint conversion amplitude of the process. Furthermore we show, 
using a simple scaling argument, that recently proposed cascaded FC protocols that overcome the 
aforementioned limitations act as “attenuators” of the TOCs. This observation allows us to argue 
that very similar cascaded architectures can be used to attenuate TOCs in photon generation via 
spontaneous parametric down-conversion. Finally, by using the Magnus expansion, we argue that 
the TOCs, which are usually considered detrimental for FC efficiency, can also be used to increase 
the efficiency of conversion in partially mode selective FC. 


Nonlinear photonic materials provide some of the most 
advanced platforms for manipulating the frequency and 
spectral profile of photons. This manipulation is typi¬ 
cally achieved using a process known as frequency con¬ 
version (fc) m 12 . In the version of FC we consider, 
two photons (one of them typically coming from a bright 
classical field) are fused into another photon with a 
higher energy by the process of sum frequency genera¬ 
tion. FC has several important applications, including 
photon detection^, and the establishment of compati¬ 
bility between sources and quantum memories [4] . FC 
can also be used to modify in a controlled manner the 
properties of weak signals or single photons, a useful com¬ 
ponent of several quantum information processing proto¬ 
cols that harness the infinite dimensional Hilbert space 
structure of the frequency degree of freedom of photons 
[SHE!- A set of orthogonal frequency amplitude functions 
for the photon(s) (henceforth referred as “modes”) pro¬ 
vides a natural basis in which to encode information in 
this Hilbert space [5], and FC provides a natural way to 
do controlled operations in this Hilbert space [B]. For all 
applications of FC it is important to have conversion effi¬ 
ciency near unity, and for controlled operations it is also 
important to have mode selectivity. In the limit of very 
short crystals, or equivalently very long pulses (effectively 
CW fields) it has been shown that 100% FC is achievable 
mm- Nevertheless in these limits most mode selectivity 
is lost. In this letter, we examine limitations to highly 
efficient mode selective FC. We show below that these 
limitations are due to time-ordering corrections (TOCs) 
that appear because the interaction picture Hamiltonian 
that describes the \2 (or X 3 ) interaction between the dif¬ 
ferent fields does not commute with itself at different 
times. Our study allows us to separate very cleanly the 
“ideal” operation of an FC device from the “undesirable” 
effects of time ordering. Thus, we can explicitly write the 
operation of an FC gate, in a language very close to the 
one used in quantum information, as a unitary operation 
ZYfc = exp(f2i + UI 2 + CI 3 + ...) that consists of a desired 
generator Hi and TOCs IF, H 3 ,... that modify the op¬ 
eration of the gate. This separation is obtained using the 
Magnus expansion (ME) jlOj. Using the ME allows us to 


understand the scaling of the TOCs as a function of the 
energy of the classical pump pulse. Based on this model, 
we propose a new scheme to achieve highly efficient (and 
partially mode selective) FC, harnessing the TOCs that 
until now have been undesirable; this model also allows 
us to explain, using simple scaling arguments, why the 
double pass scheme introduced by Reddy et al. [8] suc¬ 
ceeds in achieving high efficiency FC, and to extend these 
ideas to the case of photon generation using spontaneous 
parametric down-conversion or four wave mixing. 

We begin with the Hamiltonian describing a second- 
order nonlinear optical FC process in the unde¬ 
pleted pump approximation, occurring in a quasi one- 
dimensional structure of length L and characterized by a 
nonlinear coefficient X 2 [10] ■ 

Hi(t)= — fhe j du} p duj a dujb(e 1 ^ t ^(Ak(uj a ,ujb,ujp)L/2) 

a(w p )d(w a )S t ( w&)+H.c.), (1) 

where refers to a pump frequency, and the shape of the 
classical pump function a(oj p ) is taken to be a Gaussian, 

a(oj p ) = re~ T u r/y/Tr, Suj p = lo p — lo p , (2) 

where Co p identifies the center pump frequency; r iden¬ 
tifies the duration of the pump pulse, and we specify 
its energy by Uq. In the FC process frequency com¬ 
ponents oj a can be destroyed and frequency components 
cob can be created; associated with this are bosonic de¬ 
struction and creation operators d(uj a ) and b' (uif,) re¬ 
spectively. We assume that the phase-matching func¬ 
tion (PMF) <F(a:), where its argument Ak(uj a , U)b, ui p ) = 
kbiuib ) — k a (cj a ) — kpiojp) ± 27t/A (including the quantity 
27 t/A only if periodic poling for quasi phase matching 
is performed on the device), and the ki{uii) indicate the 
dispersion relations of the modes involved, restricts the 
destroyed and created photons to be in nonoverlapping 
frequency regions, and/or of different mode profiles or 
polarizations; then we can take [a(w a ), ¥(u>b)] = 0. We 
reference the frequencies to a and uJb to center frequencies 
w Q and Cdb respectively, defined so that energy and mo¬ 
mentum conservation are exactly satisfied for the center 


2 


frequencies, 

Cj b — u) a — w P = 0 and A k(io a ,uj b , u> p ) = 0. (3) 

Finally, A = co b — LJ a — uj p , and 


£ = 2L\2\ 


V2U 0 TTU b U} a 


^/Tr(47r) 3 eoAc 3 n a (tOa)n b (LOb)n c (uj c )7 


( 4 ) 


is a dimensionless constant that characterizes the 
strength of the interaction; the are the indices 

of refraction at the central frequencies c*>», and A is the 
effective overlap area of the spatial fields. 

Ignoring TOCs (as it is implicitly done in 1), the uni¬ 
tary evolution operator connecting states at t —> —oo to 
states at t —> oo is 

Ur = e hl = e~^ dtHl{t) (5a) 

_ g— f dui a duii,[Ji(i^a,uib)a(ui a )b^ (cj(,)+H.c.) (5b) 


where the quantity Jr is, to this level of approximation, 
the joint conversion amplitude (JCA) and is given by 

Ji(w a , u b ) = -ea(uj b - w a )$(Afc(w a ,w b , u b - u a )L/2). 


To understand how this operator transforms a given in¬ 
put state, let us first introduce the Schmidt decomposi¬ 
tion of the function J\ 

-Ji(w 0 ,w b ) = Y -^^ fc ?( w a)^(wfc), rg(e) = erg, (6) 
g 

where the Schmidt functions kg(ui) (Igfu)) are assumed 
to form a complete and orthonormal set, and the posi¬ 
tive quantities rg(e) = erg are the Schmidt numbers of 
the function Jr- In Eq. § we have explicitly used the 
fact that £ is only a multiplicative constant in this ap¬ 
proximation of the JCA; accordingly, the Schmidt num¬ 
bers are linear functions of this quantity and the Schmidt 
functions are independent of it. Because of the orthonor¬ 
mality of the functions kg(uS) and lg(ui), we can write 

yi x _ giEe re(s)(AoBl+H.c.) _ (gj giro(e)(AoBj+H.c.) , ^ 

g 

Ag = J duj a k* g (u} a )d(u} a ), Bg = J duj b ll(u; b )b(uj b ), (8) 


where we introduced broadband Schmidt operators Ag 
and Bg for fields a and b. These operators transform in 
the expected way under U\\ for example, 

U\AgU\ = cos (r g (e))Ag - ism(rg(e))Bg. (9) 

Let us consider how the unitary in Eq. ( |5a| ) transforms 
a single photon of the form, 



(10a) 

(10b) 



— Cl + 
i sin(j’i) 


+C2 + 

i sin(r 2 ) 



FIG. 1. Schematic representation of a frequency converter. 
The device has a complete set of input Schmidt modes (repre¬ 
sented in orange); an arbitrary input single photon can thus 
be represented in this basis with expansion coefficients eg. 
As explained earlier (ft, iTTj, the frequency converter acts as 
a beam splitter of reflectivity sin(j’g) in each of the Schmidt 
modes connecting input modes to output modes represented 
in blue. Note that the mode profiles can in general have com¬ 
plex amplitudes and thus we use continuous/dashed lines for 
their real/imaginary parts. 


The amplitude g is assumed to be normalized according 
to f <hj\g(u>)\ 2 = 1, and in the last equation we used the 
completeness of the set {kg(oj a )}. We can now use Eq. 
@ and Eq. ( |10b ) to st udy the effect of the unitary on 
the single photon state (10a), as sketched in Fig. [l] 


Ui\l g{uia )) - Y eg (cos(r e (£))Aj, +isin(r e (£))SjJ |vac) 
g 

( 11 ) 

In particular, if a single photon with frequency pro¬ 
file kg' (ui a ) is initially prepared, it will be upconverted 
to a photon with profile lg'(u> b ) with probability ampli¬ 
tude isin(£7-0'). One very interesting case of the last 
equation is when the JCA is separable Ji(u! a ,u)b) = 
£rofco(w a )Zg(wb) and thus has only one nonzero Schmidt 
number in Eq. (|6|. Then one can selectively upcon- 
vert only Schmidt function 0 and leave the rest un¬ 
changed; this is precisely the ideal operation of a quan¬ 
tum pulse gate (QPG) as introduced by Eckstein et al. 
[B]. However, the U\ in Eq. (5a) is not the correct uni¬ 
tary operator since very generally [Hj(t),Hi(tf)] ^ 0, 
and thus Eq. ( |5a| ) should be premultiplied by the time¬ 
ordering operator T [32], as has been mentioned in this 
context mm but not worked out in the needed detail. 
One very special an important case where Eq. (5a) does 
hold is when the PMF is flat $(A/c(w 0 , w b , w p )) ~ 1 since, 
as shown by Donohue et al. @], [Hi(t), Hi(t')] = 0 in that 
case. 

Recently we have shown that the ME provides a partic¬ 
ularly appealing way of approximating the time evolution 
operator for time dependent quadratic Hamiltonians m 
like the one in Eq. ([l]). Using the ME, the correct unitary 
time evolution operator is|10j 


Ufc = Te 


dtHj(t) _ ^ni+f22+^3 + " 


( 12 ) 


—s/— 

= C G 


In the last equation, the odd order terms are beam- 
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splitter-like operators 


„ 27 T f - 

^2n+l = — du a dLd b (j2n+l(Ua,Ub)a{u] a )b\u! b ) + H.C.) 

(13) 

and the even order terms are rotation-like operators 

^2n = — /'dw c dw(,G , 2 rl (w c , w(,)c t (w c )c(w'). (14) 
1 _ J 


c=a,b 1 


The functions J‘An+\■ G‘ 2 n can be reduced to integrals of 
the PMF and pump functions [DJI- 

One can break the RHS of Eq. (12 1 into a unitary 


involving only beam-splitter terms and one involving only 
rotation terms Ufc = ^ra^rb^bs where Uv, and U-rb, 
are generated by operators such as Eq. (jl4[) and Ubs 


is generated by operators such as Eq. (131. The lowest 


order example of this factorization is obtained by using 
the Baker-Campbell-Hausdorff formula 

Ufc = (?JraZJrb) (zJbs) ~ e^ 2 e^ 1+ ^ 3+ 5 , (15) 

where 

[O^O 2 ] J dv a du b (K 3 (uj a ,v b )a(u; a )b f (v b ) - H.c.) 
K 3 (u} a ,U} b ) _ f j, j _ G%((jJ,U a )Jl(uJ,U b ) 


The factorization Eq. (15) is useful because we can write 


the unitary for FC as an entangling operation for fields 
a and b postmultiplied by local unitaries. These local 
unitaries cannot alter the entanglement between a and b , 
and because they involve rotation operators of the form 
given in Eq. (14), they cannot change the total number of 


JCA a complicated function of to a and to b , and certainly 
introduce an imaginary component (see Eq. (|l 6 |); thus 
the kg(ui a ;£) can become complicated as well. 

Another way of achieving highly efficient FC is by elim¬ 
inating the TOCs. To this end note that the n th Magnus 
term f2 n scales with e n . Now let us imagine that instead 
of sending our single photon once through an FC device 
together with a pump pulse of energy XJq, the effect of 
which is characterized by U^c (e), we send it sequentially 
through N copies of the original FC device, each with the 
same pump pulse shape but with pump energy U 3 /N 2 
(equivalently scaling e by 1 /N). Then the time evolution 
operator is (Ufg(s/N)) N = e^ 1+ ^ 2 / JV+ ^ 3 / Ar " + -; thus by 
having more than one FC device the TOCs can be re¬ 
duced. This observation neatly explains the results of 
Reddy et aZ.[5], where highly efficient-mode selective FC 
is achieved using cascaded non-linear crystals. One com¬ 
plication with this scheme is the fact that it requires mul¬ 
tiple FC devices, which typically will imply stabilizing 
and controlling more complex interferometers. Note that 
the argument just presented makes no assumptions about 
the nature of the time dependent Hamiltonian that gen¬ 
erated the Magnus terms; thus, one can also use cascaded 
nonlinear devices to attenuate the TOCs in photon gen¬ 
eration using spontaneous parametric down-conversion 
(SPDC) as well. As noted in 10; the TOCs associated 
with SPDC have a very strong resemblance to the ones 
present in FC, and in particular the Magnus terms in 
both processes have exactly the same scaling with respect 
to the energy of the pump pulse. 

Let us now consider the possibility of achieving near 
unit efficiency FC by harnessing the TOCs that until now 
have been undesirable. That is, can we turn a bug into 
a feature? A full analysis would of course have to be 
done including the e 5 and possibly higher terms in Eq. 

(16), to either take their effects into account or determine 


photons in fields a and b. Thus the unitary Ubs contains 
all the information necessary to calculate the probability 
of upconversion. The JCA associated with Ubs is now 


J — J\ + J 3 + iK 3 + . 


(16) 


that they are negligible for the parameter space being 
considered. This involves only straightforward algebra 
and integration, but we defer a complete study to a later 
communication. Here we simply restrict ourselves to Eq. 
(16) and show how complete conversion can be achieved 


We can introduce the Schmidt decomposition of J just 
as was done in Eq. ©■ except that now, because of the 
TOCs, the Schmidt functions depend on e m and the 
Schmidt numbers are nonlinear functions of e; thus we 
must replace 


kg(io a ) —► kg(w a \E), lg(u> b ) —» lg{ui b ',e), 
Ag —> Ag(e), Bg —> Bg(e), 


(17a) 

(17b) 


in Eq. (|7). Note that Eqs ( 7|10a|ll ) will still hold after 
the rules (17]) are applied. In principle 100% efficiency FC 
could be achieved were r,;(c) = it /2 for some i and s and 
the single photon A](e)\ vac) = / duj a kg(u} a ', £)d^(o; 0 )|vac) 
sent to the frequency converter. Of course, this could 
be hard to attain experimentally, since TOCs make the 


within this approximation, and what is the physics as¬ 
sociated with this possibility. We begin by considering 
for simplicity that the JCA is separable and Gaussian 
if TOCs are ignored. This can be achieved by engineer¬ 
ing a Gaussian PMF fT4] . Gaussian functions can also be 
used as a simple approximation to the more common sine 
function one would encounter for a uniform nonlinearity 
in over a region of length L. One can match the FWHM 
of e~ lx to the FWHM of sinc(a;) by setting 7 ss 0.193. 
Assuming a Gaussian PMF, 4>(a;) = exp(— 7 a; 2 ) , we ex¬ 
pand the phase mismatch around the central frequencies 
of the fields involved to get 


&(Ak(u} a , uj b , uj p )L/ 2) = e 


_ — (sbSuJb — s a Suj a — SpSujp 


(18) 


where s t = y / 7 yL/(2v i ), ^ = dk '^ 1^=^ is the inverse 
group velocity of the i th held, and we have assumed that 
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group velocity dispersion is negligible. Under these ap¬ 
proximations, ignoring TOCs the JCA is 

Ji (u> a ,u h ) = -£T ( 19 ) 

— T H"~ (Sp Sa)(Sp ^5)5 Ma ,6 — (^p ^a,b) • 

The necessary and sufficient condition for a separable 
Gaussian J| is /i = 0. This will give an ideal QPG, 


— gi / dad^(r 0 e/(a)/(^)a(a)S t (/3)+H.c.) 


7-0 


\Z 2 ttt 
yJHaUb 1 


a = HaSva, /3 = /r b <5a; h , /(a) 



for the photon |l/( a )) = Mj|vac) = f daf(a)a^(a) |vac), 
mapping it to the photon |l/(/ 3 )) = Sj|vac) = 

f d/3f(/3)b^(/3)\vac). In Eq. (20) we have used the fact 
that the Schmidt functions of 


( |19| ) with /n = 0 are 
k 0 (uj a ) = y/jj^f*(u) and l 0 (u b ) = y/jlUfiP). To achieve 
/t = 0 and an ideal QPG one can envision two strategies. 
In the first a fixed pump pulse of duration r is assumed 
and we seek a material where the product of the relative 
inverse group velocities of the fields satisfy the constraint 
0 = p 2 = t 2 + (s p — s a )(s p — Sb)- A necessary condition 
for this to happen is that s a < s p < s b or s b < s p < s a . 
If /i = 0 is so achieved, Ji is separable and one can in 
principle upconvert photons of temporal duration /i a > r 
centered around u7 a to photons with temporal duration 
fx b > t centered around uJ b . A second strategy is to 
consider a given material with parameters s a ,s b and s p . 
Then /i = 0 i s achieved by choosi ng a pump pulse of du¬ 
ration r = yj—{sp — s a )(s p — s b ) > 0, which will convert 
photons of duration /i a to photons of duration /if,. Note 
that in both strategies the duration of the photons that 
can be transduced is always larger than the one of the 
classical pump pulse (see Eq. ©). 

Finally, let us mention that mode selectivity is lost ei¬ 
ther for a flat PMF ( s a ,s b ,s p <C r) or for a very long 
pulse (r s a ,s b ,s p ), for fi a , fj, b , fj, ss r and the Schmidt 
number S = [lalibi is very large. Nonethe¬ 
less in both cases it is possible to achieve high conversion 
efficiency [31 SO- 

As soon as one starts to approach an FC efficiency 
near unity, one needs to include the TOCs. The uni¬ 
tary connecting inputs and outputs is given by Eq. 
(15). The efficiency of the conversion is now solely gov¬ 
erned by the Wbs term in Eq. (15), and the associ¬ 
ated JCA in Eq. (16). The amplitudes in Eq. (jT 6 |) 
can be related to the terms of the TOCs for SPDC 
m via Ji(Suja L Suj b ) = Ji(-5uja,6u) b ), J 3 ( 6 iv a , 6 w,b) = 
-J 3 (5u) a ,-5uj b ), K 3 ( 6 uj a , 6 u> b ) = K 3 (-5(v a ,SuJb), where 
the unbarred quantities refer to the TOCs for SPDC. 
Note that the TOCs will typically generate more Schmidt 
functions; even if J\ is separable, J in Eq. (16) will not 
be for e ^ 1. For the very simple Gaussian PMF, and as¬ 
suming /t = 0, it turns out that the Schmid t fu nctions of 
the time ordered corrected JCA, J in Eq. (16), are only 
functions of foe. This quantity is the Schmidt number of 




FIG. 2. In A we plot the probability of upconversion as a 
function of the interaction strength foe, for a single photon 
in the Gaussian wave packet f(a) including and excluding 
TOCs. In B we plot the Schmidt numbers of J including 
(dashed lines) and excluding (full line) TOCs. In the inset 
we plot the overlap |ce | 2 between the time ordered corrected 
Schmidt functions and the Gaussian profile of the incoming 
single photon. 


the Schmidt function of interest in the e< 1 limit (see 
Eq. ©)• Also note that fo parametrizes the departure 
of the separable J\ from a circle. For a perfectly round 
JCA it takes the maximum value fo = 7 r. 

Now let us consider what would happen if a single pho¬ 
ton with fixed Gaussian profile Aj|vac) were sent to the 
FC device. We look at the probability of upconversion 
as a function of foe, and plot this probability, includ¬ 
ing and excluding TOCs, in Fig. [2}A. If TOCs could be 
eliminated (for instance following the approach of Reddy 
et al. we would get p b = sin 2 (foe); instead, we 

get a more complicated curve. The upconversion curve 
attains a maximum near the ideal foe ~ ir /2 but only 
reaches ~ 80%, a value consistent with experiments m 
The reason for this is easy to understand if we look at 
Fig. [2]B, where we plot the Schmidt numbers and the 
overlap of the Schmidt functions with the single photon 
profile /(a). We see that up to foe < 1.47 ~ 7 t/ 2 , the 
Schmidt numbers are not significantly modified by the 
TOCs. Yet the Schmidt functions are significantly mod¬ 
ified at this point, and thus do not completely overlap 
with the shape of the incoming photon. This can be 
seen in Fig. where we plot the overlap of the Gaus¬ 
sian wave packet with the two Schmidt functions. As 
we increase fo£ further, the probability of upconversion 
oscillates faster; this is just a reflection of the fact that 
past fg£ = 7 t /2 the Schmidt numbers are modified by the 
TOCs and grow much faster than linearly, as seen in Fig. 
@B. This feature could be used as direct experimental 
evidence of the contribution of the TOCs to the JCA in 
FC. Our results suggest that in some FC experiments an 
increased pump power would eventually lead to enhanced 
upconversion rates, but only after an initial drop in the 
upconversion probability. The final, most interesting fea¬ 
ture of Fig. [2] A is that the probability reaches unity for 
a value near r 3 e ~ 2.79. This happens because the two 
non-zero Schmidt numbers of J become odd multiples of 
7 t /2 at the same time. This allows the two Schmidt func¬ 
tions to cooperate and attain 100% efficiency. The extra 
Schmidt mode is, as mentioned before, generated by the 
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TOCs; thus one can think of this enhanced upconversion 
probability as an interference effect between the first and 
third order Magnus terms. Note that this upconversion 
process will be partially mode selective, since it only in¬ 
volves two Schmidt modes while leaving the other modes 
untouched. 

In this letter we have studied the effects of the time or¬ 
dering corrections (TOCs) on frequency conversion (FC) 
using the Magnus expansion. We have considered con¬ 
version using a xi material — quasi-phase-matching can 


easily included and, with correspondences presented 
earlier [.10], this can be generalized to the use of \ 3 mate¬ 
rials as well. We have shown how the TOCs modify the 
ideal operation of a quantum pulse gate and how they 
can be used to achieve near unity FC using the extra 
Schmidt functions that are generated by them. Finally, 
we provided a simple scaling argument that explains the 
results of Reddy et al. [5], in which TOCs are attenuated 
by using cascaded FC devices. 
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